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THREE-DIMENSIONAL PLASMA ACCELERATION THROUGH AXISYMMETRIC 
DIVERGING MAGNETIC FIELDS BASED ON DIPOLE 
MOMENT APPROXIMATION* 
byH. G. Kosmahl 
Lewis Research Center 

SUMMARY 

The acceleration of a tenuous plasma whose initial energy is stored in the magnetic 
dipole moments of the electrons immersed in a diverging magnetic field was investigated 
theoretically. Adiabatic conservation of the dipole moments and axisymmetric steady- 
state fields were assumed, and three-dimensional trajectories of the plasma and the in- 
duced space-charge fields were computed. The computations were carried out for atomic 
weights of 40 and 130. The particle trajectories cross the field lines, and the plasma 
emerges as a well-defined beam. From the amount of divergence in the beam, a propul- 
sive efficiency is defined. This efficiency is primarily a function of the radius and the 
length of the solenoid and the position beyond the center of the magnetic field at which the 
plasma of a given radial extension is initially injected. Conditions for a favorable design 
of a propulsion device are shown, and the mechanism of the acceleration process is ex- 
plained. Results agree well with experiments performed on a microwave -powered plasma 
accelerator working under conditions that approximate the theoretical assumptions. A 
significant conclusion of these studies is that trajectories are independent of the ionic 
mass and the magnitude of the dipole moment, that is, of the energy of the particles. 


INTRODUCTION 

The subject of this analytical study is a three-dimensional computation of electron- 
ion trajectories in a tenuous collisionless plasma. These computations are closely re- 
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The material in this report is also published in J. Appl. Phys. , no. 



lated to conditions existing in the experi- 
mental, microwave-powered accelerator 
described in reference 1. In the compu- 
tations, a number of electrons, together 
with an equal number of ions, are in- 
jected with a constant velocity u Q along 
the field lines of a cylindrical solenoid 
(fig. 1). At the plane of injection, the 
B field is still closely parallel to the 
axis of symmetry. In the model, all in- 
jected electrons are assumed to possess 
an adiabatically invariant magnetic dipole moment p parallel to B which is caused by 

the orbiting motion of the electrons around the B flux lines. Other important assump- 

1 2 2 

tions are that total initial plasma energy per particle is equal to - to£r£ • m = M 0 B res > 

6-function distribution of the dipole moments around ju = p Q , injection at a single plane 
z = z Q , single ionization of the species, absence of collisions, and zero initial angular 



'-Solenoid 

Figure 1. - Schematic diagram of accelerator. 


momentum. (All symbols are defined in appendix A. ) 

For a given set of input conditions (ft Q and (ft Q , the trajectory of a charged particle 
in steady-state electromagnetic fields is described completely by the equation 


m(ft = ±e(E + (ft x B) where (ft and (ft are the radius vector and the absolute velocity of the 
particle itself, respectively. Because of the complex nature of the motion, composed of 
the superposition of the "spiraling" (corresponding to p), the translating (r and z), and 
the azimuthal ( 6 ) components of the electron, it is easier to follow the method first sug- 
gested by Alfven (refs. 2 and 3): each particle is assigned a "substitute" trajectory de- 
scribed by the r, 0 , and z motion of its guiding center while all the "spiraling" energy 
is considered as potential energy juB. Note that the 6 motion of the guiding center must 
not be confused with the motion that produces p. 

As Hellwig has shown (ref. 4), p is an adiabatic constant of motion at least to the 
second order in the derivatives of B in spatially and temporarily changing magnetic 
fields, if, taken over one cyclotron orbit, the change AB is small compared with B; 
that is, AB « B. The validity of this assumption is examined in appendix C. 

The actual motion of an electron is obviously more complicated than that of a substi- 
tute particle because of instabilities, oscillations, and random motion. Because space- 
charge neutrality must be preserved to a very high degree over regions equal to or larger 
than some few Debye shielding lengths X^, a substitute electron will be, on the average, 
in the vicinity of each ion. In a way, the ionic inertia mass M may be considered 
tied to the dipole moment p of the substitute electron by way of space-charge forces. 
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Since any relative motion of electrons and ions in either the r or the z direction 
produces a charge separation, only the 8 motion is unimpeded. Thus, the macroscopic 
ion and the substitute electron velocities and trajectories must remain substantially equal. 
Essentially, these conditions require a macroscopic space-charge neutrality and the dis- 
appearance of the current components j = j = 0, except j 0 , which does not equal zero. 
Since electrons can acquire, in general, any velocity u < c with small energy expendi- 
tures, the expelling force on the magnetic dipole (n ■ V)B must be closely opposed by the 
induced space charge -force eE. This latter force serves, in turn, to expel (accelerate) 
the ions. 

The ejection of a plasma out of magnetic fields has also been analyzed by Kilpatrick, 
et al. (ref. 5), who treated the plasma as a conducting continuum in thermal equilibrium. 
The result of Kilpatrick's studies is that the plasma follows the B lines and cannot free 
itself from the field unless a very small angle of divergence is maintained. The results 
of the present investigation show that, while the trajectories are strongly affected by the 
B field, even a collisionless plasma can free itself from the flux lines. 

Although the problem of a three-dimensional acceleration of a plasma through axi- 
symmetric diverging magnetic fields is presented in this study as a general model, nu- 
merical computations were carried out only for cases closely associated with the 
microwave-powered accelerator (ref. 1), which might be of practical importance to 
space propulsion; that is, B vQO » 3000 gauss, co « co , « 8. 35 gigahertz, 
n(z Q ) « 10 - 5X10 X ^ particles per cubic centimeter, and u^ « 5X10 meters per second, 

corresponding to a specific impulse of 5000 seconds. 


TRAJECTORY EQUATIONS 

Trajectory equations are obtained by solving simultaneously the momentum equations, 
the energy equations, the continuity equations, and the Poisson equation. 


Momentum Equations 

The momentum equations will be derived from a Lagrangian formulation of the prob- 
lem. Cylindrical coordinates and azimuthal symmetry are assumed throughout. The 
generalized Lagrangian L* of the guiding center can be written for electrons and ions, 
respectively, in the form 



( 1 ) 
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( 2 ) 


L* = A Mu? 
1 2 1 



with the notations 


u = a r r + a 0 r0 + a z z (3) 

<P e = -V + (4a) 

e e 

(p i = -V (4b) 

For the electrodeless, steady-state model to be investigated, the following conditions are 
applied: 


j r = j z = 0; E 0 = -grad 0 V = 0; -1 = 0 (5) 

0t 

A = a ft A; B = curl A = a JB_ + (6) 

C7 r r z z 

In equation (6), any magnetic fields induced by the dipoles have been neglected as small 
compared with the external field produced by A. 

The momentum equations are obtained from 


£ 

dt 




(7) 


where q and q are the coordinates and the rate of charge of coordinates, respectively, 
and and are the corresponding gradient operators. From the assumption that 
H is adiabatically invariant and — = — = 0 equations (7) and (1) yield the momentum 
equation for electrons: 


u. 



/m - v y _ _ 

[-y + » e x(7xA)j=i 

*V-(-7- 2 ) B + ' 1 e XB 


( 8 ) 


Because of the assumption of the adiabatic invariance of p, all the derivatives of /i are 
considered small compared with the derivatives of B and, consequently, are neglected 
in equation (8) and the following equations. From equations (7) and (2), the momentum 
equation for the ions is 
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( 9 ) 


£ = - [-V-V + U. X (v X A)] = -1 (-V V + £ X 5) 
1 M L q 1 J M v q 1 ' 


The term 
of F are 



represents a force, which shall be denoted as 


F. 


The components 


3B_ 


3B, 




F„ = — - + 0 + u 

r r 3r z 3z 


Fq = 0 + 0 + 0 


3B_ 3B_ 

F z = M r — - + 0 + \i 

z r 3r z 3z 


( 10 ) 


with 


M = -M 


B(r, z) 
° |B(r,z)| 


u = constant > 0 , 
o J 


(ID 


The induced dipole moment is antiparallel to the magnetic field and is thus negative. 

With the help of equation (10), the r, 0, and z components of equations (8) and (9) 
become, for the electrons. 


F 

r = r'9 2 - ± (E r + B z r0 e ) + -J±L (12a) 

m m 

- - PO = - - (B r z - V* (12b) 

r dt V / m 

■m 

z = - — ( E z - B r r0 e ) + (12c) 


and for the ions, 
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r = rflf + ± <E r + B z re,) 

(13a) 

- i(r 2 8i)= +-< E z i - B i-) 
r dt V M z z 

(13b) 

z = + — (E - B re.) 
M z r 1 

(13c) 

Next, in integrating equations (12b) and (13b) use is made of Busch's theorem in axisym- 
metric fields (ref. 6). For electrons and ions, respectively, one obtains: 

m(r 2 0) e = +e[i//(r, z) - ^ 0 (r Q , z Q )] 

(14a) 

and 


M(r 2 &)j = -e[^/(r, z) - ^ 0 (r Q , z Q )] 

(14b) 

where 


i//(r, z) = B z (r, z)r dr 

(15a) 

^o (r o’ z o> = ^ °B(r,z 0 ).rdr 

(15b) 

2* 2 m 

The initial angular momentum mr Q 0 o or Mr o 0 o of the guiding centers in equations (14) 
has been set equal to zero, in accordance with the assumption that the dipoles and the ions 
are injected with zero angular velocity. Equations (14) are used to eliminate 6 from 


equations (12) and (13) in the computations. 


Energy Equations 

The energy equations are obtained by multiplying equations (8) and (9) by u g or u^ 
respectively: 
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( 16 ) 


IA(u)2 = + ±Sv.J.p:2 = ± A (y - iLIU - -l ^ 

2 dt e m dt m dt m dt \ e / m dt 

1 i. (u.) 2 = - - dV - e 

2 dt 1 M dt M dt 

Equations (16) and (17) are easily integrated to 

j [u 2 (r, z, e Q ) - u 2 J = p- [B Q (r 0 , z Q ) - B(r, z)] - e[v(r, z) - V^r, z)] (18) 

and 

M [u 2 (r, z, e Q ) - u 2 ] = e[v(r Q , z Q ) - z)] (19) 

For z — °°, B w and grad B M as well as grad V approach zero. The plasma is 
assumed to have a reference potential = 0 at z — 


Continuity Equation 


The continuity equation is developed with the help of figure 2. As a consequence of 
zero net current, j z = j r = 0 and macroscopic space-charge neutrality p net ~ 0, u e (r, z) 
must be always close to or identical with u.(r, z). Thus, u(r, z) w u. (r,z) « u.(r, z). In 

X v 1 

figure 2 a current density ±P' 0 U 0 (the signs indicating ions or electrons, respectively) is 
injected perpendicular to the plane z = z Q between r Q and r Q + dr. The current dI Q 

through this annular ring of area 2nr o dr Q is 



d! 0 = ±2ffp 0 u 0 r 0 dr 0 (20) 


To calculate dl, the current between r and 
r + dr perpendicular to DC, it must be noted 
from figure 2 that 


tan a = 


r’ 


DC = dr • sin a = 


dr 

Vl + (r') 2 


where 
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CALCULATION OF BEAM DIVERGENCE EFFICIENCY 

If particle trajectories are assumed known throughout the entire exhaust region in- 
cluding infinity, a » ’divergence' * efficiency tj d , which is defined as the ratio of the di- 
rected kinetic energy of the stream at infinity to the total kinetic energy, can be calcu- 
lated from 


8 


(24) 



N 



1 


2 2 
z + r 


w 


1_ 

N 


N 

V-is- 

Z-/ 1 + (r*) 2 
1 


where N represents the number of particles injected per second. That there is no azi- 
o. 

muthal energy Mr flj at infinity should be noted. 

The function f introduced in equation (24) is a weighting function representing the 
number density of particles injected initially between r Q and r 0 +dr Q per unit time such 
that 

r D / 2 

IX V.V‘ < 25 > 


Two cases are examined here. First, if the propellant injection density n Q at z = z Q 

were constant across the waveguide, there would be 27 m o u Q r o dr Q particles injected per 

second between r and r + dr . Thus, f„, = n^u^. 

o o o ’woo 

Since r'(r 0 ) is known along the trajectory and is continuous across the injection 
plane, the summation in equation (24) may be replaced by an integral and the following is 
obtained for rj: 


%26 " 



(26) 


The second case assumes that the propellant is injected preferentially along the axis 
r«0 in such a manner that n Q = — , where a is a proportionality constant. The effi- 
ciency becomes r o 


%27 " 


2i r 
N 


1 


f r <r 1A 


r <r° 


1 + 


^o 


2ir 

N 


a 



4 


dr o 


(27) 
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The actual propellant injection density in the quoted experiments lies somewhere between 
the two cases. 


SOLUTION OF EQUATIONS 


The preceding equations (12), (13), (14), (15), (18), (19), (22), and (23), constitute 
in a mathematical sense a complete and uniquely defined set of expressions for a given 
external magnetic field. Although this set of equations is solvable in principle, the gen- 
eral behavior and the order of magnitude of the solutions must be determined a priori in 
order to make the computational program tractable. This problem is dealt with in detail 
in appendix B. 

The range of parameters and numerical values chosen for the computations corre- 
sponds to that for the accelerator of reference 1. The driving microwave frequency f is 
8. 35 gigahertz. Consequently, a peak magnetic field B > 3100 gauss is necessary for 
achieving the resonance condition « u> c . Peak values of B were chosen in these 
computations to permit the resonance to occur at z Q = 2. 5 and z Q = 5.0 centimeters. 

The static magnetic field B was calculated from the expression for the vector poten- 
tial A of a single layer solenoid of radius R and length L wound only in the 6 direc- 
tion. The equation for A q is (ref. 7) 


and 


A e 


A e 



I 1 (kR)K 1 (kr)cos kz sin — 
1 1 2 


K^kRjI^krJcos kz sin — 
1 1 2 


dk 

k 


dk 

k 


for r > R 


for r < R 


(28a) 


(28b) 


Thus, 


B(r, z) = V X Aq 

The constant C is determined from the condition that at the origin 

r = z = 0, B(0, 0) = B q 
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Most experiments were carried out with argon and xenon, with mass flow rates be- 
tween 0. 5 and 2.0 milligrams per second and power levels between 500 and 2000 watts. 

From these numbers and the atomic masses M = 6.68X10 kilogram and 
-26 *** 

Mxenon =21. 95X10 kilogram the injection rate in particles per second N lies between 

18 19 • 17 

10 and 10 atoms per second and an energy per particle P/N lies between 5X10” 

and 2X10"*^ joule. 

The value u can then be calculated from u » P / — by assuming a 6 function 

O O 

res 12 

distribution of p and no energy losses. The variation of u is between 2X10 and 

12 0 
60X10” ampere - square centimeter. Because of the accent on low specific impulse, 

-12 -12 

the range selected for computations was approximately 2x10 to 8X10 ampere - 
square centimeter. 

The plasma "diameter" at the injection plane is identical with the diameter D of the 

circular waveguide in the TE^ mode for an X -band frequency range; that is, D = 2. 8 

centimeters. If it is assumed temporarily that the injected plasma fills the waveguide 

cross section evenly and that it is injected with a sonic velocity u at room temperature. 

• 2 ^ 

The propellant density n Q can be obtained from n Q « 4N/ (7rD • u Q ) and amounts to a few 

12 2 2 
times 10 per cubic centimeter for u = 1.7X10 to 3X10 meters per second. The 

u 12 

limit for n Q (ci>p ^ oj^) is approximately n Q * 5X10 . 

The computational step in computing the trajectories was 0. 1 centimeter for r and 
0. 5 centimeter for z. 

The entire computational loop used in solving the equations consisted of the following 
major steps: 

(1) B -field, equations (28a), (28b), and (6), and the dipole force field F (eq. (10)) 
were mapped. 

(2) Flux integrals (15a) and (15b) were computed. 

(3) First approximations to trajectories from the condition if; = constant were found. 

(4) Laplace and Poisson equation (23) were solved on axis. 

(5) The potential V was mapped off axis. 

(6) A more accurate trajectory was computed. 

(7) A more accurate V-field was computed. 

DISCUSSION OF RESULTS 

The results of the computations are presented in figures 3 to 5. Figure 3(a) shows 
the trajectories for argon and xenon injected at z Q = 2. 5 centimeters, and figure 3(b) 
shows the same trajectories for an injection at z Q = 5.0 centimeters. Figure 4 shows 
fields and potentials for various values of )J. Q - All figures were computed for a 
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Radial position, r, cm 


(a) z 0 - 2.5 centimeters. (b) z Q - 5.0 centimeters. 

Figure 3. - Trajectories of seventh (r n - 0. 7 cm) and fourteenth (r n - 1.4 cm) argon and xenon particles for three solenoid geometries. 




Potential V, V 


\ 

B X \ V 

B \ 


v-r = 1.5 
X.- r = 0 


13 14 

13 7 


\ 

\\ \ 

\\ 




(a) \i Q = 2. 08xl0" 12 ampere-square centimeter; z 0 * 2. 5 centimeters. 
E z on axis and E z at r - 1. 5 centimeters. 


- E z ( r - 0) 


0 20 40 60 

Axial position, z, cm 

(b) V on axis; Mq = 2. 08xl0' 12 ampere-square centimeter; z 0 : 
2. 5 centimeters. 


Figured - Magnetic field B, induced voltage V, induced radial 
electric field E r , and induced axial electric field E z against dis- 
tance z for three solenoid geometries. E r at r = 1.5 centimeters. 
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(d) p 0 * 7. 82xl0' 12 ampere-square centimeter; 
Zq ■ 5 centimeters. 


Figure 4. - Concluded. 


magnetic field of peak value B q « 4000 
gauss, but for three different dimen- 
sions of the field producing solenoid: 

R = 7, L = 7; R = 13, L = 14; and 
R = 13, L = 7 centimeters. 

The following conclusions can be de- 
duced from a review of the results: 

(1) For equal R, L, and z Q and the 

same initial plasma injection radius, the 

macroscopic trajectories are independent 

of the ionic mass M and the magnitude of 

the dipole moment /i (fig. 3). The range 

° - 12 

of p. was approximately 2X10 to ap- 
-12 

proximately 8x10 ampere - square 
centimeter. 

(2) The trajectories depend on R, L, 
and z Q , and their slope r’ is increasing 
with decreasing R and L for a fixed 
plasma injection diameter. Comparing 
figure 3(a) with figure 3(b) shows, in addi- 
tion, that the slope r* also increases 
with decreasing z Q . 

(3) The induced voltage V is always 
near -proportional to the magnetic induc- 
tion B (fig. 4). 

(4) The induced electric field E z is 
almost independent of the radial position r 
(fig. 4) and peaks in the neighborhood of 

z « L/2 (figs. 3(a) to (c)). The maximum 
of E, and of E„ vary with the gradient 
of B and, thus, vary inversely with the 
solenoid dimensions R and L. 

(5) For a fixed magnet geometry and 
field the induced potentials and electric 
fields are independent of the ionic mass M 
and proportional to p, which in turn is 
determined only by the ratio P/N, in ac- 
cord with the assumptions of a monoener- 
getic plasma. The difference in the calcu- 


14 


.4 


.020i 



0 10 20 30 40 50 60 70 


Axial position, z, cm 

Figure 5. - Slope r' and derivative of slope r" of fourteenth 
particle (r 0 -1.4 cm) against distance z for solenoid R = 13 
centimeters, L = 7 centimeters; zq = 2.5 centimeters. 

lated potentials for the same ju results from local differences in the value of B in the 
injection plane and from a different location of z Q of the plane itself (figs. 3(b) to (d)). 

(6) The plasma exhaust velocity is approximately equal to u^ » y2/i 0 B res /M and 
depends only on the ratio ju/M, but not, however, on \i and M individually, as long as 
M » m. 

(7) Figure 5 is a plot of r' and r" for the outermost particle as a function of z 
for one selected case demonstrating the convergence of the trajectory solutions. As r' 
converges asymptomatically to its final value, r M approaches zero. In figures 3(a) 
and (b), the trajectories are approaching straight lines at axial distances z ^ 20 to 40 
centimeters. 

(8) The discontinuity in E in figures 4(a) to (c) at z = z^ results from the discontinu- 
ous introduction of the magnetic dipoles and is equivalent to a sudden space-charge force 
acting at z Q . In a real accelerator, such a discontinuity is "smoothed” out and E is 
continuous. 

(9) The azimuthal velocity of the ions is less than 3 centimeters per second 

* m 5 6 

everywhere, and that of the electrons r'9 = — r 9. is on the order of 10° to 10 centi- 

C ^ i 

m 
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TABLE I. - EFFICIENCIES 


Initial axial 
injection 
position, 

z o> 

cm 

Ratio of 
radius to 
length, 
R/L, 
cm 

Efficiency, 77 

Equation (27) 

Equation (28) 

2.5 

13/14 

0.974 

0.964 


13/7 

.971 

.960 


7/7 

.920 

.880 

5.0 

13/14 




13/7 

0.975 

0.964 


7/7 

.943 

.90 



meters per second. Both velocities ap- 
proach zero at large values of z in 
accordance with equation (14). 

(10) The efficiencies, as defined in 

equations (27) and (28), are tabulated in 

table I for z =2.5 and 5. 0 centimeters, 
o 

three R, L values for the solenoid, 
and a waveguide diameter D = 2. 8 cen- 
timeters. The loss of directed beam 
energy resulting from the divergence 
of the plasma exhaust is not critical as 
long as 2R/D s 5*. A noticeable effect 
of z Q on i) also occurs. The im- 
provement of 7] with increasing z Q is 
the result of much less entrapped ini- 
tial flux for a given solenoid, so that 
less divergence occurs. 


COMPARISON WITH 
EXPERIMENTAL DATA 

The available experimental data 
which were obtained with the acceler- 
ator pictured in figure 6 permit the 
following comparisons with theory: di- 
vergence angle of the exhaust and its 
independence of the ionic mass M, 
slope of the induced space -charge po- 
tential V, and the magnitude of V. 

Figure 7 shows typical photographs 

of the visible plasma exhaust produced on a photographic plate by argon and krypton pro- 
pellants. The solenoid geometry was R = 13 centimeters and L = 7 centimeters. (While 
the luminous beam edge itself is very clear, its photograph is not because of fluctuations 
during the exposure time. ) Of particular interest is the divergence angle of the outermost 
trajectory (r = 1.4 cm, beam edge), which is drawn in figure 7 as the angle y. In good 
agreement with the exit angle y in figure 3, which is approximately 18°, r Q = 1.4 centi- 
meters, R = 13 centimeters, L = 7 centimeters, y is found to be approximately 19° to 


Figure 6. - Electron cyclotron resonance accelerator; 8. 35 gigahertz. 
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(a) Argon. Gas flow rate, 0.93 milligram per second. 



(b) Krypton. Gas flow rate, 1.48 milligrams per second. 

Figure 7. - Exhaust stream. Coil current, 210 amperes; power, 
1 kilowatt radio frequency. 



Accelerator exit to calorimeter entrance distance, cm 


Figure 8. - Dependence of calorimeter power on calorimeter 
position. Propellant, xenon; calorimeter, 10 inch diameter 
by 20 inch length. 


Curve Mass flow, Accelerator exit to 



mg/sec 

array distance, 
cm 

1 

1. 07 

72 

2 

.87 

112 

3 

.76 

132 

4 

.97 

192 


1 



I I I i I i I i I 

0 8 16 24 32 

Radial position, cm 

Figure 9. - Exhaust stream ion current density profiles at 
several axial positions as determined by diode tube array. 
Propellant, argon; radiofrequency power, 684 watts. 
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20° in figure 7 and is quite independent of the ionic mass. Additional evidence (ref. 1) 
from calorimetric sampling of the exhaust (fig. 8) and ion density measurements (fig. 9) 
shows that the luminous beam is near-representative of the actual exhaust stream. 

From figure 4, it was theoretically concluded that V is proportional to B. Figures 
10(a) and (b) are experimental plots of the axial magnetic field and the induced plasma po- 
tential for argon and xenon, respectively. Over the measured region, V is near- 
proportional to B. In figure 10(b), the measured potential flattens out at z > 25 centi- 
meters instead of following B. This deviation, however, is the result of an excessive 



12.5 17.5 22.5 32.5 42.5 77.8 12.5 17.5 22.5 32.5 42.5 77.8 

Distance from window, cm 


(a) Argon. (b) Xenon. 

Figure 10. - Magnetic field and plasma potential profiles against distance. Radiofrequency power, 1 kilowatt. (Numbers from 
operating conditions listed in Table VII, ref. 1). 
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background pressure (ref. 1) due to diffusion and charge -exchange processes which tend 
to short out V. 

The magnitude of the induced space-charge potentials should, according to figures 
4(a) to (c), range from approximately 150 to approximately 500 volts at axial positions 
between z = 15 centimeters and z = 40 to 50 centimeters. From figures 10(a) and (b), 
on the other hand, only a potential drop of about 50 volts can be deduced. This large dis- 
agreement is not yet fully understood. Processes that could affect the calculations in- 
clude background pressure effects, diffusion processes, ionization and excitations, insta- 
bilities, oscillations, a distribution in ju Q , and all direct energy losses. Furthermore, 
the driving microwave is known to be imperfectly polarized and its E field is not con- 
stant across the waveguide. Also, other modes than the "desirable” TE^ must appear 
in the vicinity of the waveguide edge, where a transition from a bounded TE ^ wave into 
a free space wave occurs. Finally, ionization and energizing of the dipoles do not occur 
suddenly at a single injection plane z Q , but rather take place over a finite distance on 
the order of a few centimeters. 

CONCLUSIONS 

The following conclusions can be drawn from the results of a study of three- 
dimensional plasma acceleration through axisymmetric diverging magnetic fields: 

1. The particle trajectories do cross the flux line, and come out in a well-defined, 
slightly diverging beam . 

2. The trajectories depend on solenoid geometry and position of injection, and their 
slope increases with increasing plasma cross section relative to solenoid dimensions. 

3. To keep thrust losses resulting from beam divergence small, the plasma diameter 
should not exceed one-fifth of the solenoid diameter, and the injection plane should be 
moved out axially toward the magnet face. 

4. For the same solenoid geometry, the same plane of injection relative to the cen- 
ter of the magnetic field, and the same plasma diameter, the trajectories do not depend 
on the ionic mass and the dipole moment of the plasma electrons. The induced electric 
fields are a function of the electron injection energy only. 

5. The induced space -charge potential is proportional to the product of the magnetic 
field times the dipole moment and is only weakly dependent on the radial position. 

6. To obtain physically sensible solutions, the initial enclosed flux must be larger 
than the flux at an exit plane. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 13, 1966, 

120-26-03-00-22. 
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APPENDIX A 


SYMBOLS 


a r’ a z’ a 0 


B 

®res 

C 

c 

D 

E 

e 

F 

f w 


W 


Ki(k r ) 

k 

L 

L* 

M 

m 


N 


magnetic vector potential 
unity vectors 

curl A, magnetic field/ 

T = 10 8 G 

magnetic field for electron 
cyclotron resonance 

computation constant 

speed of light, 3xl0 8 m/sec 

defined in fig. 2 

gradient of V, electric field, 
V/cm 

1. 6X10“ 48 coulomb 
magnetic force on dipole, N 
weighting function 
current, A 

modified Bessel function of 
first kind and first order 

9 

current density, A/cm 

modified Bessel function of 
second kind and first order 

Fourier wave number 

length of solenoid, cm 

generalized Lagrangian 

mass of ion 

rest mass of electron, 

9. lixio" 28 g 

ion (atom) injection rate, 
particles/ sec 


n 

P 

P r 

q, q 

R 

(ft 

r c 

r, z,0 
s 

t 

u 

V 

V 


a 

6 

€ o 


particle density, cm 
power 

power radiated in cyclotron 
radiation 

coordinates, rate of change of 
coordinates 

radius of solenoid, cm 

radius vector of particle 

electron cyclotron radius 

cylindrical coordinates 

path element 

time 

a r r + a 0 re + a z z - velocity, 
cm/sec 
potential, V 

transverse velocity of electron 
in cyclotron orbit 

r (? e - azimuthal electron ve- 
locity 

injection constant 
Dirac 6 function 
8. 86X10“ 14 F/cm 
divergence efficiency 

magnetic dipole moment, 

(A) (cm 2 ) 

space-charge density, 
coulomb/cm 

2tt/ u> c 
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<p 

generalized potential, V 

* 

f B z (r, z)r dr, 

•'O 


magnetic flux/2 n 

U) 

frequency 

w c 

cyclotron frequency, e/mB 

"p 

4 

plasma frequency, 5.6X10 

"rf 

radian frequency, cps 

Subscripts: 

e 

electron 

eff 

effective 

i 

ion 


net neutrality 

o initial 

res resonant 

r, z, 0 cylindrical coordinates 

00 final 

Superscripts: 

total time derivative, 

3-3 

— + u 

at 3s 

' partial spatial derivative, 3/3z 

M second spatial derivative 
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APPENDIX B 


APPROXIMATIONS FOR INITIATION OF COMPUTER PROGRAMS 


In this section, a number of approximate expressions will be developed and order of 

magnitude estimates will be established for the initiation of a computational program on 

the IBM 7094, in order to make it tractable. 

The character of the force F on the electrons (eq. (10)) should be considered first. 

aB z 9B z 

In the vicinity of the axis F_ « p • and is positive since both u and are 

z z 9z 9B z dz 

negative. The predominant term in F_ is u , which is negative over the initial 

9z 

most important range of acceleration. As defined in the equations the electrically in- 
duced field E„ and E„ is considered positive if it accelerates the ions in the positive z 
z r 

and r direction. Subtracting equation (13a) from equation (12a) and neglecting cor- 
responding terms divided by M as compared with those divided by n yields 


r e - r i " rS e - eE r 


B r0 + — « 0 
z e 

m 


(Bl) 


.. . o 

because r = r” • z + r’z must be almost identical for ions and substitute electrons so 
that the left side of equation (Bl) is approximately zero. The symbol y will be substi- 
tuted for r 6 e and treated for the moment as the only unknown in equation (Bl). The 
quadratic equation 


has two roots 





(B2) 


(B3) 


For r on the order of a few centimeters the term 



6 r 2 

r - B z 

m 


is on the order of 10 


6 


be- 
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cause E r and F r /e are both less than 10 volts per centimeter. Thus, very accurately 



(B4a) 



(B4b) 


The first term in equation (B4a) is (for B on the order of kilogauss) substantially larger 

than the second term and the kinetic energy - my^ can assume any arbitrarily large 

2 + 

values, which are determined only by the amount of B. This situation cannot be true, in 

general, and equation (B4a) must therefore be rejected. Notice also that the total energy 

of the plasma in our model cannot exceed m~B„_ 0 per particle. Equation (B4b), on the 

o res 

other hand, admits any arbitrarily small value for y and represents, therefore, the only 
acceptable solution. 

As mentioned before, F r < 0 in the region where B and therefore E are large, the 

last term in equation (B2) is therefore positive, and the first term y^ is positive but small 

because of our choice of equation (A4b). The term r — B y = — B r^6 must therefore 
2 . mm 

be negative, or r *6 < 0, which is possible only if, in equation (14a), 

iMr, z) < ^ 0 (r Q , z Q ) (B5) 

The inequality equation (A5) is of fundamental importance , for it shows that the plasma 

must follow a tube of a slightly decreasing flux in order to obtain physically compatible 

solutions. If (F r /e) - E r is looked upon as being an effective electric field E r ef £, then 

equation (B4b) may be regarded as the 6 component of E^xB/B^ drift. At larger 

F 

values of z, F r is no longer negative; however, E r remains positive and | E r | > — , 

e 

so that the last term in equation (B2) is less than 0 and the compatibility conclusion made 
earlier is still valid. Thus, the plasma trajectories do not strictly follow magnetic flux 
lines, and because they do not, the plasma can free itself from magnetic fields. 

Addition of equation (14a) and (14b) leads to 

mr^ e + Mr^0j = 0 (B6) 

which means that the total angular momentum ctf the plasma is zero if it was zero at the 


23 



injection plane. This latter assumption of zero initial angular momentum is very reason- 
able for the case considered here as well as for a number of other cases with realistic 
injection schemes. 

Equation (B6) determines, because r very closely equals r., the angular momen- 

Q , v 1 

turn of the ion species after r 0 g = ry has been calculated from equation (B4b). The 

calculations of E„ and E„ are discussed next. On the axis E = B =0. Since the 
z r r r 10 

space -charge densities p Q and p. in equation (23) must be equal to within 10" or less 

in order that the induced voltages be around a few kilovolts or less, the ion and electron 

velocities must almost be equal everywhere. Thus, from equaling equations (12c) and 

(13c) the following is obtained: 

E = (B7) 

z e e B 9z e 3z 

on the axis, and E_ » F/e elsewhere as long as B„ < B„. For the latter condition, 
z z r z 

3B B 3B 

F » p = , and if a similar argument is applied to equations (12a) and (13a), 

3z | B | 3z 

the following may be written as a fair approximation, which is better than an order of 
magnitude: 


E, 


B, 


E. 


B. 


(B8) 


which assures the continuity of E r — 0 on the axis, r = 0. Notice that equation (B8) is 
only a first-step approximation in the computation process and therefore is not identical 
with the final computed values. 

Integration of E_ along the axis gives V(z) and z(z). Integration of 

n 

/ E(z = constant, r dr) gives V(r, z = constant) with the initial value V(0, z = constant) 
J r=0 

used as the constants of integration. Mapping of the complete potential field V(r, z) de- 
termines the ion velocity u^(r, z) = u^ + z^|l + (r’)^J + r^0? from equation (19), with r’ 

yet to be determined. This value of r f can be estimated with the help of equation (B4b). 
Since the radius of the trajectories is between 0 and about 0. 1 meter, the order of mag- 
nitude of r 0 is 10 square meters per second. According to equation (14a), r 6 Q is 
equal to e/m times the flux difference ^ - i l/ Q . Again, for B « 0. 3 tesla and 

r < 0. 1 meter, — i// » 10^ and — ip « r^0 + — \j/ « 10^ + 10 » 10^ is obtained. Thus, 
m m m ° 

the substitute electron must follow a trajectory that is close to a tube of constant flux. 

The latter is easily computed as long as the induced magnetic field is small compared 
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4 

with that externally applied, and an initial approximation of r' is that which results 
from constant flux trajectories. 



APPENDIX C 


DISCUSSION OF BASIC ASSUMPTIONS 

In this section, the validity of some basic assumptions made in this study will be 
checked. 

The first assumption concerns the adiabatic conservation of /i. That AB incurred 
during one cyclotron orbit is much smaller than the local value of B will be proved at a 
position where B is small and thus at the least favorable condition. Select B = 500 
gauss; then from B field mapping is found 

M = 50 G/cm 
3z 

AB = — x Az = H u x T = 50 X 4. 5 x 10 6 x 1 x 10“ 9 = 0. 15 G, 

3z 9z c 3 

6 

if stream velocity u is approximately 4. 5x10 centimeters per second and 

w = — B « 9X10 9 sec -1 
c m 


and 


t = 2L » 2. x 10 -9 
C w c 3 


Therefore, 


AB _ 1.50X10 _ 0 3 x 10 -3 

B 500 


which is negligibly small. Thus, ju is conserved down to even much smaller values of 
B, at which the acceleration process is completed. 

Next, it will be proved that the power radiated in cyclotron radiation is small for the 
conditions assumed here. The cyclotron radiated power is (ref. 8) 


P . 2 e4 Z 4 B 2v i . ! e « XB 3 y 
r 3 „ 3„2 3 __ „3 3 


4ire Q c m‘ 


7re 0 c m 


(Cl) 
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because v 2 = and Z = 1 where Z is the charge number; for /! « 8X10" 12 ampere - 
m 

1 O 

square centimeter and B = 3000 gauss, P„ is 6x10" watt per electron. With 
* 18 r 

N = 3. 65X10 atoms per second, the energy input per particle is 1000 watts per 
3.65X10*8 « 3x10 *8 joules. The acceleration time is approximately 

50 cm * 50 cm « 20X10" 6 sec 

it max 2. 5x10® cm/sec 
2 

Even if the electrons radiated at maximum power (they do not because B decays), the 
radiated energy would be less than 6xi0”*8x20xl0"8 * 1.2X10" 22 joules « 3X10"*® 
joules. Equation (Cl) is valid for Vj^ « c, which is correct for the range of parameters 
selected for this computation. 
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